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Flatness-based control of a single qubit gate 

Paulo Sergio Pereira da Silva and Pierre Rouchon 



Abstract 

This work considers the open-loop control problem of steering a two level quantum system from 
an initial to a final condition. The model of this system evolves on the state space 3£ = 51/(2), having 
two inputs that correspond to the complex amplitude of a resonant laser field. A symmetry preserving 
flat output is constructed using a fully geometric construction and quaternion computations. Simulation 
results of this flatness-based open-loop control are provided. 



Index Terms 

Quantum control, nonlinear systems, geometric control, flatness, qubit gate. 

I. Introduction 

Take a single qubit, i.e. a two level quantum system. Denote by COq its transition frequency. 
Assume that it is controlled via a resonant laser field v G M: 

v = wexp( — i(OQt) +K*exp(zcooO (1) 

where u = u\ + 1U2 G C, (wi,«2) £ ^ 2 > is its complex amplitude. In general, the frequency COq is 
large and the time variation of u is slow: \u\ <C (Oq\u\. In the interaction frame, after the rotating 
wave approximation and up to some scaling (see e.g., [1]), the Hamiltonian reads u\0\ +U2O2 
where 0\ and 02 are the first two Pauli matrices (see appendix). The gate generation problem 
then reads: take a transition time T > such that (OqT <C 1 and a goal matrix U G SU (2); find 
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a smooth laser impulsion [0,T] 3 t i— > G C with w(0) = u(T) = such that the solution 
[0, T] 3 t U(t) G SU(2) of the initial value problem 

ijU(t) = (u 1 (t)a 1 +u 2 (t)a 2 ) U(t), U(0)=I 2 (2) 

reaches U at time T, i.e., U(T) = U. This motion planning problem admits a well known 
elementary solutionQ. It relies on the fact that U = exp(— iJO\) exp(— i[5a 2 ) exp(— i(XO\), for 
all U G SU(2), for convenient (a,/3,y) G M 3 (see, e.g., [2]). An obvious steering control u(t) is 
decomposed into three elementary and successive pulses: for the first (resp. third) pulse, u 2 = 
and u\ is such that its integral over the pulse interval equals a (resp. 7); for the second pulse, 
u\ = and the integral of u 2 is /3. 

Here, we propose another solution where u\ and u 2 vary simultaneously, i.e., the steering 
control u(t) is contained in a single pulse. Our solution does not rely on optimal control 
techniques (see for instance [3] and the references therein) and is explicit. It does not rely 
on numerical resolution scheme. It provides controls that can be chosen to be C a or C°° function 
of t. As far as we know, such explicit solution is new and could be of some interest for reducing 
the transition time T while still respecting the rotating wave approximation. Our approach is 
based on the fact that the system dynamics is differentially flat [4]. The flat output constructed 
in this paper has a clear geometrical interpretation. 

In section HH theorem \T\ shows, using a quaternion description of ©, that this system is 
flat. We propose a coordinate free definition of the flat-output that lives in the homogenous 
space SU(2)/ exp(iRoi). This geometric construction preserves invariance with respect to right 
translations. In the sense of [5], the flat output is compatible with right translations. The proposed 
construction can be seen as the analogue of the geometric construction based on the Frenet 
formula for the car system, where the steering angle is directly related to the curvature of the 
path followed by the flat-output curve [6]. In section [Till we show how to use such geometric 
flatness parameterization to solve analytically the motion planning problem corresponding to 
an arbitrary quantum gate. Simulations illustrate theorem [2] and the interest of such explicit 
open-loop steering control. In section [TV] some conclusions are briefly stated. Some material 
has been deferred to the appendix . In part A one finds the basics properties of Pauli matrices 

'The so-called ZYZ quantum logic gate. 



February 2, 2008 



DRAFT 



3 



and their associated quaternions as well as the correspondence between SU (2) and quaternions 
of length one. In part B one finds a proof of the fact that the motion planning algorithm has no 
singularities. 

II. A SYMMETRY PRESERVING FLAT OUTPUT 

The dynamics © reads in quaternion notation (see appendix A) 

^-q= {u\e\ +u 2 e 2 )q (3) 
at 

where q E Hi is a quaternion of length one and where («i, W2) £ ^ 2 is the control relative to the 
modulation of a coherent laser field (u\ + 1U2 is the complex field amplitude). This system is a 
driftless system on the Lie Group Hi. It is controllable (see, e.g., [7]). Moreover, this control 
system is invariant with respect to right translations in the sense of [5], [8]: 

• the group G = Hi acts on the state space = Hi via right multiplication g : q i—> qg 
where q 6 Hi. 

• the dynamics is G-invariant: if 1 1— > (q(t),u\(t),U2(t)) is a solution of © then 1 1— > (q(t)g, 
u\{t), U2{t)) is also a solution of © for any g e G. 

The controllability structure of this system is in fact of a very special kind. Around any point 
q e Hi, © can be seen in local coordinates as a driftless controllable system with 3 states^ and 
2 controls. Thus, as known since [9] (see also [10]), such system is differentially flat and the 
flat output function can be chosen to depend only on the state. More precisely, the flat output 
for the controllable system ^x = uifi(x) + «2./2(X) with dim(jc) = 3 is obtained by the rectifying 
coordinates of any vector field f(x) = (X\(x)f\(x) + which is a linear combination of 

the two control vector fields f\ and fi (CL\, OL2 are any scalar functions of x). 

We propose here a coordinate free and symmetry preserving construction of the flat output 
via the previous procedure. Thus we are looking for a flat output map h : Hi 1— > where <3f 
is the output space, a compact manifold of dimension 2, and G-compatible in the sense of [5]. 
This means that the output map h must satisfy the following constraint: there exists an action of 
G = Hi on the flat output space W described by the transformation group p g : y 1— > p g (y) such 

2 Take e.g., the exponential map: (xi,X2,*3) £ R 3 i-» exp(;qei +-*2£2 + x 3 e 3) 1 ? mat ma P s a neighborhood of EE R 3 to a 
neighborhood of q in Hi . 
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that p g (h(q)) = h(qg) for any q E Hi. The following construction will be based on the control 
vector field associated to u\, and hence to e\. 

Denote by K = {exp(^>ei)}^ e [ ,2^] me one dimensional subgroup of Hi generated by e\. 
We can consider the action of K on Hi via left multiplication: to any k E K, we have the 
diffeomorphism Hi 3 q \— > kq E Hi. Two elements of Hi, q and p, belong to the same orbit if 
and only if there exists k E K such that kq = p. Denote by W the set of the orbits. This set is 
a compact manifold of dimension 2 and the output function h is the map that associates to any 
q, the orbit h(q) to which q belongs. This map is a smooth submersion, and & is called an 
homogenous space (see, e.g., [11]). If q and p belong to the same orbit, qg and pg also belong 
to the same orbit for any g E Hi. Therefore, this output map is G-compatible in the sense of [5]. 

Assume that y(t) is a curve on W = H\/K, at least of class C 2 . Since the map h : Hi — > W 
is a submersion, in adequate local coordinates one has h{x\,X2 : x^) = (xi,X2). Assume, without 
loss of generality, that the open neighborhood of definition of h is rectangular and contains 
(0,0,0). Define locally the smooth map g : U C W — > V C Hi, where U, V are open sets and 
g(x\,X2) = (xi,^2,0). Note that g is smooth, and Y(t) = g(y(t)) is such that h(Y(t)) =y(t). 
Then, locally, there exist smooth maps g^ and g^ such that Y(t) — g^\y(t),y(t)) and Y(t) = 
gW(y(t),y(t),y(t)). 

Let us show now that the map h defines a flat output. This means that the inverse of system © 
with output y = h(q) has no dynamics^. Thus we have to consider the following implicit system 

^q = (u\e\ +u 2 e 2 )q, y = h(q) 

where 1 1— > y(t) is a known function of time and where the quaternion q(t) G Hi and the control 
(u\(t),U2(t)) are the unknown quantities. 

The problem is how to manipulate h, since only a geometric construction for h is available. 
Knowing the function 1 1— > y(t) means that we have at our disposal a smooth function 1 1— > Y (t) E 
Hi such that y(t) = h(Y(t)). Hence, to have y(t) = h(q(t)) means that q and Y belongs to the 
same orbit for each time t. Therefore, there exists k(t) = exp(0(?)ei) in K such that q = kY. 

3 This is equivalent to say that the state q and the input u = (u\ , 112) can be written respectively as q = (y,y,y, . . . ,y^ a ') and 
u = 33(y,y,y, . . . ,y^) for convenient smooth maps si and S3. 
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Since k{t) = q(t)Y*(t), then k(t) is smooth. Thus, we have 

—q=t—k )Y + k—Y. 
dt \dt J dt 

But jjk = Q)e\k where ft) = ^0. Using ©, we get the following equation k^Y = {{u\ — (0)e\ + 
u 2 e 2 )kY ', that is 

k ( ~r Y ) Y * k * = ("l - a) ) e \ +u 2 e 2 . 



dt 

This quaternion equation gives in fact k as a function of (jjjY) Y*. Left and right multiplication 
by ^3 yields 

e-ik I — Y J y*£*e 3 = (u\ —CO)e\-\- u 2 e 2 



dt 

since £3^3 = for i = 1,2. Hence, we have the following relation (without the controls and 
ft)): 

e3 ^yWe 3 = ^W (4) 



dt 



Assume that 

d 



— Y )Y = CQiei + G>2,e2 + a>3e3 (5) 
dt 



where the ft) ; 's are known smooth real functions of time. Thus, we get 

k \^-Y) Y*k* = (0\e\ + k 2 (o> 2 e 2 + G*se-$) 



dt 

since e\k* = k*e\, kk* = 1 and = kei for / = 2,3. Therefore, © reads: 

k 4 (c02e 2 + ah,e 3 ) = (co 2 e 2 - e 3 e 3 ) 

since e^k 2 = (£*) 2 <?3 and k~ l = k*. 

Right multiplication by e 2 yields the following algebraic equation defining k 

fc 4 (ft)2 + ft) 3 ei) = (ft>2 - ft^ei). 

Since k = cos0 + sin0ei, we have the following equation for the angle 

(cos40 + sin40 ei)(a> 2 + ft> 3 ei) = {(O2 — ft^ei) 

2 

which is equivalent to exp(40i) = A? where z = ft>2 — ftty is a known complex number. Thus, 
there are four distinct possibilities for k: 



fe = ±exp^eij, k = ±e\ exp Q^ei J 



(6) 
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where is the argument of 6O2 — 6O3Z. The controls u\ and U2 associated to one of these four 
trajectories are obtained by 

e^k—YY*k*ej ) = (u\ — co)e\ -\-u2e2 
where 2(0 = is given via simple algebraic formulae based on 6O2, 6O3, ^6)2 and ^603: 

6O3 £ CO2 -CO2^C0 3 
60 = — — 

2(cq| + (0 3 2 ) 

For the two branches k = ±exp (f ei) we get 

CD 3 |ft)2-ft)2|0>3 

Mi = 6O1 H — = ^ 

2(60 2 + 60 2 ) 



and for the two other ones k — ±e\ exp (f ei) we get 



u\ = (0\ H — 2 — 2- — 

2(601 + 601) 



w 2 = -^/©l + col 

where just the sign of j/2 is changed. All the previous computations are valid when 6O2 — 6032 ^ 0, 
i.e., when 4v 7^ 0: (60 2 + 60 2 ) does not depends on Y(t) such that h(Y(t)) = y{t); it depends only 
on y{t) and vanishes if, and only if, 4- t y{t) = 0. To summarize, we have proved the following 
result: 

Theorem 1: Take 7 > and an arbitrary C 2 curve [0,7] 9 ? 1— »y(f) on ^ such that 4fy{t) 7^ 
for any t G [0,7]. For any smooth curve f 1— > F(f) G Hi such that /z(y(f)) = set z = 
602(0 _ G>3(t)i 7^ for all t G [0,7] where (^Y) Y* = (i>\e\ + 602<?2 + Then there exists 

a smooth function [0,7] 3 1 1— > 6(t) G R such that exp(6h) = A and any smooth solution 1 1— > 
(q{t) 1 u\(t) 1 U2{t)) of © satisfying h{q(t)) =y(t) for all ? G [0,7] is one of the four following 
trajectories indexed by n G {0, 1,2,3}: 

I (t) = (e l fcxp( ( ^ C1 W) 



Ll>V 

Ml = 60l + 



t0 3 iC02 - 602^603 (?) 



2(601 + 601) 
u 2 = (-l) n yjG$ + a% 
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Recall that some Y(t) such that y(t) =h(Y(t)) is locally given by Y(t) = g(y(t)), and fur- 
thermore Y(t) = g^\y,y) and Y(t) = g^ (y,y,y), where g, g^ and g^ are smooth maps. In 
particular, the last theorem proves that y = h(q) is a flat output. 

The flat output y — h(q) is obtained with e\ playing a specific role. In fact one can see that any 
map h-q : Hi i — ► Mi/K^ (rj £ S 1 ) corresponding to the subgroup = exp(M(cosr]ei +sinr]e2)) 
also defines a flat output. It just corresponds to a rotation by the angle Tj of {q\,q 2 ) and («i,M2)- 
If we set 

e\ = cosT]ei +sinr]e2, ^2 = — sinr/ei +cosT]e2 

the imaginary quaternions (^1,^2^3) and (S\,e2,e^) satisfy exactly the same commutation rela- 
tions. Thus, if 1 1— > g(f) is a solution of © with the control (mi(?),M2(0) men 

? ^>qo(t) + (cos r] ^x(f) - sin T]q 2 (t))ei + (sin r] ^1 (?) +cos7]^ 2 (0) e 2 + ^3(0 e 3 

is also a solution of © with the control 

u\ = cosT]«i(?) — sinTji^^)) "2 = sinr]Mi(?) +cosr]«2(0- 

This symmetry and the fact that, as stated in theorem [H h — ho is a flat output, implies directly 
that h-q is also a flat-output. The family (/ijj)jj e §i is made of flat outputs all compatible versus 
right translations. 

III. Motion planning 

In this section, we will use C7]) with n = to propose an explicit solution for the motion 
planning problem stated in the introduction: for any 7 > and any final state q £ Hi, find a 
smooth control [0, 7] 3 1 1— > u(t) = (ui(t),U2(t)) £ K 2 with w(0) = w(7) = 0, such that the solution 
[0, 7] 3 1 1— > q(t) £ Hi of © starting from q(0) reaches q at time 7: i.e., q(T) = q. 

As the system is driftless, every time re-parameterization of a solution is also a solution. In 
fact, consider the equation 

= (ui(s)ei + U2(s)e 2 )q(s) 

Let g : [0,T] — > [0, 1] be an increasing diffeomorphism. Then q(s) is a solution of the previous 
equation defined on [0,1], with input (u\(s) 7 U2(s)) if and only if q(t) = q(g(t)) is a solution 
of (|3]) defined on [0,7] with input («i(?),H2(0) = Tjf ("l(?(0)>"2(?(*))- O ne concludes that, 
without loss of generality, one may always state the motion planning problem with the (virtual) 
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time s belonging to the interval [0, 1] and after that, one may "control the clock" by choosing a 
convenient bijection s = g(t). Thus, it is enough to solve the motion planning problem in the s 
scale where we can disregard the fact that the control has to vanish at the beginning and at the 

2 3 

end: it is enough to take for example q(t) = 3 (£) — 2 to get u equal to zero at t = and 
att = T, since |g(0) = f t g(T) = 0. 

In the sequel we propose a solution in the s-scale. For clarity's sake, we will remove the ~ 
when u and q are considered as function of s. The derivation in s will be denoted by ': du/ds = u' , 
dq/ds = q' , ... 

Thus, we have to find a smooth control [0, 1] 3 s \- > u(s) such that the solution of 

q'(s) = (ui(s)ei + u 2 {s)e 2 )q{s), q(0) = 1 

satisfies q{\) = q, where q is any goal state in Hi. 
We can always assume that 

q = qo + \Jq 2 l +ql(smf)ei+cosf)e 2 ) +q~3 

for some angle fj e [0,2^]. Thus, as explained at the end of last section, up to a rotation of 
angle r) of the control, we can assume that q\ = 0. More precisely, if q\ ^ 0, set rj to be 
the argument of the complex q 2 + qii- If s t— > (u\(s), u 2 (s)) steers q from g(0) = 1 to g(l) = 
<7o+ \J^ + Q 2 ^2 + G3e3, then the control 

s i— > (cosf)ui(s) +smf)u 2 (s), — sin f)u \(s) +cosf)u 2 (s)) 

steers q from g(0) = 1 to g(l) = q. 

Thus up-to a rotation of angle f) of the control, we can assume that q\ = and q 2 > 0. Thus 
we can define two angles a e]0, tt] and j8 e [— f ,f] such that 

<? = cos a + sin a (cos j8 ^2 + sin j8 63). 

If the control 5 1— > u(s) steers the system from g(0) = 1 to <7( 1) =cos a + sina(cosj8e2 + sin/3<?3), 
the same control steers the system from 

q(0) = cosA + sinA(cosj8e2 + sinj8e3) 

to 

q(l) — cos(X + d) +sin(X + Ct)(cosj5e 2 + singes). 
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This is a direct consequence of right translation invariance of © and right multiplication by 
cos A + sinA(cos/3e2 + sin/3e3). 

Take now the formulae © in the s- scale with 

Y(s) =cos(a(s)) + sin(a(s))(cos(P(s))e 2 + sm(P(s))e3 (8) 

where a(s) and j8(s) are smooth functions such that 

a(0)=l, a(l)=A + a, jS(0) = jS(l) = jS. (9) 

Set, as in theorem Q] 

Y'Y* = (Oi(s)e\ + oyi{s)e 2 + 0)3(5)^3. 
Simple computations shows that 

z= €02 — 0^1 = exp(— iP)(a' — j/3'cosocsina). 

Now we shall construct © such that q(s) = exp(0 (s)ei)Y(s), s £ [0, 1] is a trajectory of the 
system. We will assume that q(0) = Y(0) and q(l) =7(1). So we must have 0(0) = 0(1) = 0. 
Furthermore, if we can ensure that s 1— > z(s) never vanishes, and 0(0) = 0(1) =0, then the 
trajectory of © with n = will provide a steering control u. 

Let us now show in detail how to design the functions a(s) and j}(s) satisfying these 
constraints. First of all we have the initial and final constraints ©. By taking 

-f, forae[f,f]; 
f — fj otherwise; 

we always have cos a sin a far from when s = and 5=1. Thus we can impose the following 
initial and final constraints for /3': 

p'(p)=-J^L js'(i)- asin ^ 



sin A cos A sin(A + a)cos(A + a) 

and for a' 

a'(0) = a'(l) =cos/3a. 

Then a(s) and j8(*) are the polynomials of degree < 3 satisfying these initial and final constraints. 
Since a > and |j8| < \ , st— > a(s) can be a strictly increasing function on [0, 1] and a' > for 
s g]0, 1[ (see appendix B). Thus the complex number 

z = exp(— ifi) (a' — z/3'cos a sin a) 
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never vanishes for s g]0, 1 [. For s = and 5 = 1, we have 

a' — z/3'cosasina = exp(?/3)a. 

Thus z(0) = = 6c > 0. To summarize the closed path [0, 1] 3 s i— > 2(5) g C never passes through 
nor turns around 0. We satisfy the assumption of theorem \T\ in the s-scale. Moreover we can set 
z(s) = r(s)exp(id(s)) with r(s) > and 6(s) smooth functions on [0, 1] with 0(0) = 0(1) = 0. 
We avoid with such design of a (s) and fi(s) the monodromy problem associated to the resolution 
of (exp(z0)) 4 = z 2 /\z\ 2 . Finally we have proved the following result. 

Theorem 2: Take q = qo + qie\ -\-q~2e2 +<?3C3 G Hi with q 7^ 1. Chose f] G [0,27r[ such that 
gi^i +<?2^2 = \/<?i + q 2 (sinf}e{ +cosf)e2)- Define ae]0,Jt] and /3 G [—§,§] such that 



<7o + +<?2 e 2 + <?3^3 = cosa + sina(cos/3e2 + sin/3e3). 



Set A = -§ if a G [f , ^] and A = f - § otherwise. Define a(s) and jS(j)) as being the unique 
polynomial functions of degree < 3 such that (' stands for d/ds) 

a(0) = A, a(l)=A + a, a'(0) = a' (I) = acos/3 
J 8(0) = j8(1) = j8 

jg/(Q) = _ a sin/3 pl{V asin/3 



sin A cos A sin(A + a)cos(A + a) 

Define (0\(s), 6>z(s) and 6>$(s) by 

(Oi = (l-2cos 2 (a))/3' 

CO2 — 1CO3 = exp(— if})(a' — j/3'sinacoscc). 

Then CO2 and 0)3 never vanish simultaneously and the control 




I) 



-g(0 [ COST? sinr? I ( Wl+ ^M±^ 

^ » -sinT] cost] / I , / co 1 + co 2 



s=g{t) 



steers system © from q(0) = 1 to q(T) = q with 1 1— > g(t) G [0, 1] being a C increasing bijection 
between [0, 7] and [0, 1] k> 1. When in addition -^\ s = for s = and 5=1, and n = 1, . . . ,k, 
the control 1 1-> is C^ 1 with = for s = and 5 = 1. 

Figure \T\ illustrates the steering control described by theorem [2] with T = 2, go = e^, and 
£(f) = 3(t/T) 2 — 2(t/T) 3 . We see that the control is a smooth function with maxima around 
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Fig. 1. The steering control and trajectory derived from theorem [2] with T = 2, q = and 

g(t) = 3(t/T)2 — 2(f/J)3. The control magnitude is very close to an ZYZ control design with two 

separated ^ pulses. The simulation code (matlab m-file and scilab sci-file) can be downloaded from 

|http : //cas ■ ensmp . f r^rouchon/publications/PR2 007/CodeMatlabScilabQubit . zip| 



n/2, a value close to the ZYZ design based on two successive pulses: (i<i,i<2) = (0, §) for 
t G [0, 1] and {u\,U2) = (f ,0) for t G [1,2]. Thus our flatness based design yields, with the same 
transition time and control magnitude, smooth control actions. 



IV. Concluding remarks 

The results of this paper holds if the laser matches exactly the resonant frequency. If we have 
a frequency offset of A r from resonance, then this offset leads to the following drift (see, e.g., 
[1]): 



dt 



q = {u\e\ + U2e2+A r e3)q. 



It is still interesting to notice that h(q) is also a flat output. In this case, the key relation © 
becomes 
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and k = exp(^ei) is a root of the following polynomial 

£ 4 (C02<?2 + Ct>3<?3) +2k 2 A r e3 - (©2^2 - ft^) = 0. 

Then one could try to apply similar techniques for solving the motion planning problem for this 
system, although the time-scale s = £(t) cannot be considered in this case. 
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Appendix 

A - Pauli Matrices and Quaternions 
The Hermitian matrices 

/ o 1 \ / — / \ /10 

<?i = , 02 = , 03 = 

\ i o y yi J yo -\ 

are the three Pauli matrices. They satisfy of = 1, OkOj = —OjOk f° r k 7^ j, and 

OiG 2 = io 3 , o 2 a 3 = iOi, o 3 Oi=io 2 . 



February 2, 2008 



DRAFT 



13 



Any matrix U in SU (2) reads 

U = qo - q\iC\ - qii02 - ^3^3 

with (^0,^1,^2,^3) € such that q\ + q\ + q\ + q 2 = 1. By setting 

e\ = -iO\, e 2 = -i<J2, e 3 = -i<J3 

on can identify SU (2) with the set of quaternions 

q = qo + qiei+q 2 e 2 + q3e3 

of length one. This set is denoted by Hi and corresponds to quaternions q eW such that qq* = 1 
where q* = qo — q\e\ — q2^2 — q3^3 is the conjugate quaternion of q. Thus the dynamics © 
becomes © with q corresponding to U . Notice that Hi is a compact Lie group of dimension 3. 
We recall here some useful relations for k = 1,2,3, j 7^ k and 6 E: 

e| = -l, e k ej = -eje k , exp(^)e k ) = cos^) + e k sin(j) 

exp(Qe k )ej = eyexp(-0^) 

ei<?2 = e 3 , ^2^3=^1, e 3 ei=e 2 

5 - Proof that z = CO2 — 1CO3 never vanishes for s e]0, 1[ 

Since (O2 — zd>3 = exp(z/3)(ct / — 2/3' sin a cos a), it suffices to show that a' > for j g]0, 1[. 
For this, let 8 = a — a'(0) = 6c(\ — cos/3) > 0. A simple exercise shows that the polynomial 
a(s) =as 3 +bs 2 + cs + d meeting the restrictions a'(O) = oc'(l) and oc(l) — Cf(0) = a is such 
that a = -28, b = 38, c = a'(0) and d = a(0). In particular a'(s) = -68s(s - 1) + a'(0). If 
cos/3 ^ 1, then -68s(s - 1) > 0, for s e]0, 1 [. As ct'(O) > 0, then a' > for 5 e]0, 1 [. If cos/3 = 1, 
then 8 = and ot'(0) = acos/3 > 0. So a' > for s e [0, 1]. 
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